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RESEARCH SUMMARY 


The Stochastic Life Cycle Model (SLCM) simulates 
the life cycle of anadromous salmonids and is designed 
to mimic the basic mechanisms regulating populations 
of Pacific salmon, while capturing some of the intra- 
annual and interannual variation inherent in these 
populations. This model was designed for population 
viability assessments combining advanced modeling 
techniques with concepts from the field of conservation 
biology. 

While the basic structure of the SLCM is similar to 
other life-cycle models, it differs in several ways. First, 
the SLCM incorporates stochastic or probabilistic 
processes at each step in the life cycle. The binomial 
distribution is used extensively to introduce demo- 
graphic stochasticity in survival; the beta distribution 
is used to introduce environmental stochasticity. 
Because of its stochastic nature, the model’s predic- 
tions must be expressed in probabilistic terms. Multiple 
games are run using a Monte Carlo approach to 


generate probability distributions for future outcomes. 
Second, the model is designed to use inputs from 
more detailed models for specific life stages, in com- 
bination with a minimum number of empirically based — 
parameters. SLCM users can choose among alterna- 
tive models for the more contentious life stages, such 
as juvenile migration and adult harvest, incorporating 
the results of their preferred models. An ancillary 
calibration model has been developed that allows the 
SLCM to be fitted to a historical time trace of popula- 
tion estimates, constraining expectations of survival 


_ and their variances to historical levels. 


The model also allows considerable flexibility in 
describing the dynamics of juvenile production. Users 
can choose among three density-dependent relation- 
ships to describe egg-to-smolt survival, including the 
Beverton-Holt, the Ricker, and a logistic response 
function, or use empirically based conditional probabili- 
ties. A variety of scenarios involving hatchery and 
natural production are possible, ranging from natural 
production only, to a combination of hatchery and 
natural production involving supplementation of adults, 
fry, or smolts. Allocation of naturally produced and 
hatchery-produced adults among terminal harvest and 
hatchery and natural spawning follows a set of adjust- 
able rules that affords priority to natural escapement 
and hatchery broodstock needs. 

The model is written in the SAS® programming 
language, which allows the model to operate on a 
variety of computing systems and provides enhanced 
flexibility in the analysis of model output. Users of the 
model need not be proficient in SAS. Simple input 
forms and ancillary programs to analyze model results 
allow users to run the mode! with a minimum of prior 
instruction. 
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INTRODUCTION 


Models are commonly used as exploratory tools to help fishery managers 
and planners estimate the potential effects of alternative harvest, mitiga- 
tion, and enhancement strategies on target fish populations. Examples in- 
clude Goodyear (1977), Taylor (1981), Walters (1981), Jensen and Hamilton 
(1982), Peterman (1982), and MacCall and others (1983). Within the 
Columbia River Basin, the Northwest Power Planning Council relies 
heavily on the System Planning Model (McConnaha 1992) for analysis of 
mitigation and enhancement opportunities for Pacific salmon and steel- 
head. Like many other models used in fisheries science, the System Plan- 
ning Model is a deterministic model that is designed to mimic population 
dynamics. Within a deterministic model, simulated populations rigidly fol- 
low mathematical rules that govern survival through each life stage and 
from one generation to the next. Intrapopulation variability plays no role; 
interannual variation can be introduced only by changing annual inputs 
such as river flows. 

The caricature of reality portrayed in a deterministic model is sufficient 
to guide management decisions in some cases, but can be dangerously mis- 
leading in others. To understand why, first consider their conceptual un- 
derpinnings. The more commonly used models of fish population dynamics 
were derived from surplus production models developed by fisheries scien- 
tists in the 1950’s (Beverton and Holt 1957; Ricker 1954). These models 
were derived from the gross dynamics of major fish stocks observed over 
several years. The general response of these stocks to increased levels of 
fishing harvest led pioneering researchers to postulate simple models de- 
scribing the population dynamics observed in the field. Paradoxically, 
these models have found wide acceptance among fisheries scientists despite 
the models’ generally abysmal fit to empirical data (Rothschild 1986). 

At best, deterministic surplus-production models capture only the cen- 
tral tendencies of a population. When combined with models that link 
management actions to model parameters, deterministic models can pro- 
vide a normative ranking of alternative strategies if three conditions hold. 
First, the model must be valid for the population in question, that is, the 
model should faithfully mimic the population’s underlying dynamics. Sec- 
ond, the population must be relatively robust so annual deviations from the 
expected behavior do not have long-term repercussions. And third, model 
parameters must be estimated accurately. 


In a sense, we can never be certain that these conditions hold true. We 
can be more confident with populous stocks where sufficient data provide 
a convincing argument for the model’s use. Unfortunately, this is not the 
case for many stocks of concern, such as those within the Columbia River 
Basin and other coastal streams of the Pacific Northwest. 

Addressing the uncertainty inherent in fisheries management requires 
us to restructure our perspective on planning and the use of models 
(Walters 1986). Some basic modifications in model design can help us cope 
with our inability to unerringly predict future outcomes. The first step to- 
ward a richer analysis that embraces uncertainty is to use stochastic mod- 
els. Stochastic models can incorporate at least three sources of uncer- 
tainty: (1) temporal variation in population structure and environmental 
conditions, (2) intrapopulation variation among individuals, and (3) uncer- 
tainty in parameter estimates. All three sources of uncertainty can have 
important policy implications. 

Nowhere is uncertainty in planning more critical than for stocks that 
have a high probability of becoming extinct or losing significant genetic 
resources through declining populations. Chance occurrences can be cata- 
strophic for threatened or endangered species on the brink of extinction. 
Central tendencies or expected values are insufficient when dealing with 
these populations. The probabilities of catastrophic outcomes must be as- 
sessed using stochastic models that can simulate uncertainty. Determinis- 
tic models that consider only central tendencies have no place in the analy- 
sis of threatened or endangered species. 


FUNDAMENTAL PRINCIPLES 


Stochastic models are created in two basic ways. One option is to start 
with a deterministic model and recast the model parameters as random 
variables drawn from selected probability distributions. This option intro- 
duces little or no change in the basic model structure, but does require 
specifying probability distributions for each model parameter. Particular 
attention must be given to parameter estimation. The estimation proce- 
dure must explicitly include the proposed variance structure, and it should 
produce parameter estimates with meaningful statistical properties. For 
nonlinear models, the expected or mean values of model outputs such as 
population size will differ from the output of the deterministic model using 
mean parameter values. This result, which arises from Jensen’s Inequality 
Theorem, illustrates the fallacy of assuming that the results of determinis- 
tic models using point estimates will be directly comparable to those of a 
properly constructed stochastic model, or that one can convert a determin- 
istic model to a stochastic model simply by adding noise to the parameters. 

The deterministic model with random parameters is not very satisfying 
conceptually because the deterministic relationships continue to be empha- 
sized. Stochasticity in the model is essentially noise obscuring a determin- 
istic signal. At the level of human resolution, however, nature is inher- 
ently stochastic, not deterministic. Therefore, it seems appropriate to 
build this inherent stochasticity into our models from the ground up. 

A more proper way to incorporate nature’s uncertainty is to use stochas- 
tic process models. In this approach, the focus is on the state variables in 
the model, rather than the parameters. Given the state of the system at 
time ft, the likelihood or probability of all possible future states at time 
t+1 are assessed. The range of possible future states together with their 


probabilities defines a probability distribution that is the fundamental 
building block of a stochastic process model. 

A simple example can illustrate the differences among the different 
types of models. Consider the survival of a group of individuals from time 
t tot+1. Let N(t) denote the number of individuals alive at time ¢, and 
N(t+1) the number alive at time f+1. A simple deterministic relationship is 


N(i+1) = N@)p (1) 


where p = the proportion surviving through one time unit. If N is re- 
stricted to integer values, then similar restrictions must be placed on p as 
well. For example, if N(¢) = 10, then p = 0.5 is reasonable. But if N(¢) = 9, 
then p cannot be 0.5 because that will result in N(t+1) = 4.5 individuals, 
an unrealistic outcome. Simply rounding N(¢+1) to the nearest integer will 
not work. As long as p 2 0.5, the population will never die out—as t = 
0,1,2, ...00, N(¢) = 10,5,3,2,1,1,...,1. In this model, the difference between 

p = 0.5 and p = 0.4999 is the difference between immortality and certain 
extinction. 

This simple deterministic model could be made stochastic by defining a 
probability distribution for p. We would have to define at least one more 
parameter to describe the variance of p and also consider many of the same 
integer-counting problems that are associated with the deterministic 
model. While there are practical solutions for this simple, one-step ex- 
ample, the problems mount as more steps and more parameters are added. 

The more elegant and parsimonious way to address this problem is to 
view it as a stochastic binomial process. In the binomial model, each indi- 
vidual is assumed to have the same probability of surviving (p). N(f+1) is 
viewed as the sum of N(¢) independent trials with possible outcomes 0 or 1 
(0 = individual dies, 1 = individual survives). This gives rise to the familiar 
binomial probability distribution 


Prob {N(t+1) = X} = —NO! px [1 — pyro (2) 
[N(@)—X]!X! 


which is defined for 0 < X < N(¢) and has mean = N(¢)p and variance = 
N(t)p[1- p]. The binomial model is widely applicable as a model of death 
processes in ecology (Pielou 1977) and is used extensively within the 
SLCM. In the following discussion, the independent variables N and p 
used in the binomial distribution are assigned to a variety of state vari- 
ables and probabilities, but their interpretation relative to the binomial 
distribution is the same. 

Two extensions to the binomial, the multinomial and the binomial-beta 
distributions, are also used within the SLCM. The multinomial is appro- 
priate when there are more than two possible outcomes, all outcomes are 
mutually exclusive, and all outcomes have fixed probabilities of occurring 
that sum to one. For example, the allocation of subbasin escapement 
among returns to the hatchery, subbasin harvest, natural spawning, and 
natural mortality follows a multinomial. 

The binomial-beta distribution is used where the probability of survival 
(p) is not assumed to be constant, but rather is itself a random variable. 
The variation or uncertainty in p can result from three sources: (1) tempo- 
ral changes in the environment or environmental stochasticity (Gilpin and 
Soulé 1986), (2) intrapopulation variability among individuals arising from 


genotypic and phenotypic heterogeneity, or (3) uncertainty about the 
proper value of p resulting from experimental error. Press (1989) describes 
how the beta distribution arises as a model of uncertainty surrounding p in 
a binomial process, given limited information. 

In the binomial-beta distribution, a beta distribution defines the prob- 
ability distribution of p. Using the notation introduced above, the 
binomial-beta distribution is defined by 


Prob(N(t+1) =X} = —N@! __ BilX+al, IN@-X+b]) ie) 
[N@-XIX! Bia,b) 


where a and b are parameters of the beta distribution, and 6 {a,b} is the 
beta function. The moments of the binomial-beta are 


mean = WN(t)a/[a+b] 
variance = N(¢)abla+b+N(¢)] / [((a+b)2(a+b+1)]. 


Boswell and others (1979) provide a fuller discussion of the derivation of 
the binomial-beta and its use in ecology. 

While natural stochasticity and measurement error can be expressed in 
identical terms mathematically, they might have different implications 
when the results are interpreted. Thus, it may be useful to keep the two 
forms of uncertainty separate in an analysis. The SLCM has the capability 
to accommodate many such instances (see the operational notes below). 


OVERVIEW OF MODEL 


The SLCM can be thought of as a series of compartments corresponding 
to stages within the life cycle of anadromous salmonids. The total popula- 
tion is divided into two stocks, hatchery and natural. The natural stock in- 
cludes all fish spawned in the wild and hatchery-produced juveniles that 
are released as fry, regardless of the origin of their parents (fig. 1). Simi- 
larly, the hatchery fish include all fish spawned in the hatchery and re- 
leased as smolts (fig. 2). For much of their lives, hatchery-produced fish 
and naturally produced fish share a common life cycle. Transitions from 
one compartment to the next determine the model dynamics. At each tran- 
sition, random draws from the appropriate probability distribution deter- 
mine the values of the state variables (table 1). The model runs on an an- 
nual time step, but multiple transitions can occur within a single step. 
Each transition is described in more detail in the following sections. State 
variable names are given in capital letters and parameters in italics. 

The model uses a set of production, passage, and harvest parameters 
(table 2) gleaned from the various input files (described below). A set of 
control parameters determines the number of replicates (games) within 
each simulation, and the number of years per game. The user must specify 
initial numbers for each life stage. As the model runs, it creates a data 
set containing identification variables and updates state variables at the 
end of each year. For simulation purposes, a “year” begins with natural 
spawning. 

Currently, two versions of the base model are available, incorporating 
two different approaches to modeling production of juvenile migrants 
(smolts). Version 1 is more descriptive, as it tracks egg production and 
egg-to-smolt survival in a more conventional form. The user can choose 
among three different functional forms to describe density-dependence in 
early juvenile survival. Interactions between hatchery-released fry and 
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Figure 1—Flow diagram for naturally produced fish in the 
Stochastic Life-Cycle Model. Interior boxes denote 
probability distributions used in each step. 
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Figure 2—Flow diagram for hatchery-produced fish in the 
Stochastic Life-Cycle Model. Interior boxes denote 
probability distributions used in each step. 


Name 


Natural stock 


SPAWNERS 
FSPAWN 
EGGS 
PRESMOLT 
SMOLTO 
SMOLTi 
SMT_BON 
RECRUITS 
OiHARV 
OiRVH 
OiSBE 
OCHARV 
RVHARV 
SUBESP 
SUBHARV 
TOTHARV 
MAXHT 
HATTAKE 


Hatchery stock 


HSPAWNER 
HFSPAWN 
HEGGS 
HATFRY 
FRYREL 
HSMOLT 
HSMT_BON 
HRECRUIT 
HOIHARV 
HOIRVH 
HOISBE 
HOCHARV 
HRVHARV 
HSUBESP 
HSUBHARV 
HTOTHARV 
HMAXHT 
HSPAWN 
ESPWILD 


Table 1—State variables used in the Stochastic Life-Cycle Model 


Description 


Spawning adults 

Female spawning adults 

Eggs spawned in the natural environment 
Juvenile fish immediately before migration 
Subyearling juvenile migrants 

Juvenile migrants of age /+ 

Smolts surviving to below Bonneville Dam 
Fish harvested in fisheries or escaping to subbasin 
-ocean fish harvested in ocean (1<i<4) 
-ocean fish harvested in Columbia River 
ocean fish escaping to subbasin 

Total ocean harvest 

Total mainstem harvest 

Total escapement to subbasin 

Subbasin harvest 

Total harvest 

Number available to :the hatchery 

Number taken by hatchery 


Spawning adults 

Female spawning adults 

Eggs spawned in the natural environment 
Number of fry produced 

Hatchery juveniles released as subyearlings 
Hatchery juveniles released as smolts 
Smolts surviving to below Bonneville Dam 
Fish harvested in fisheries or escaping to subbasin 
ocean fish harvested in ocean (1<i<4) 
-ocean fish harvested in Columbia River 
-ocean fish escaping to subbasin 

Total ocean harvest 

Total mainstem harvest 

Total escapement to subbasin 

Subbasin harvest 

Total harvest 

Number available to the hatchery 

Number taken by hatchery 

Number escaping to natural spawning areas 


Name 


Natural production parameters 


ffem 
prspsv 
egfem 
stdegg 
jJackspn 
logtscl 
alpha 
beta 
cvegsv 
stayi 
inbsmsv 
basercv 
cvadtrv 
minesp 
maxthar 


Hatchery stock parameters 


hffem 
hprspsv 
hegfem 
hstdegg 
egg_ftry 
hfry_sv 
fry_smt 
hsmt_sv 
hjackspn 
hbaserv 
hevadtrv 
hstray 
htermhar 


Hatchery management parameters 


adtsup 
colrate 
hcolrate 
minrun 
hatrem 
hatneed 
htarget 
eggcap 
fryplant 
fryint 


Adult recovery parameters 


(h)ocnihar 
(h)ocnirvh 
(h)ocnisuv 
(h)ocnisbe 


Passage parameters 


(h)msurv 
(h)passcv 


Control variables 


ngames 
years 
smttog 


Table 2—Biological parameters and control variables used by the Stochastic Life-Cycle Model 


Description 


Probability of being female 

Prespawning survival probability 

Mean eggs per female 

Standard error of mean eggs per female 
Spawning contribution of one-ocean fish (0-1) 
Logistic response scale parameter 
Density-independent juvenile production parameter 
Density-dependent juvenile-production parameter 
CV of egg-to-presmolt survival 

Conditional probability of staying over ith winter 
Subbasin presmolt-to-smolt survival 

Fraction recovered as adults in calibration years 
CV of basercv 

Minimum escapement goal 

Maximum terminal harvest (number) 


Probability of being female 

Prespawning survival 

Mean eggs per female 

Standard error of mean eggs per female 
Egg-to-fry survival rate 

Initial subyearling survival following release 
Fry-to-smolt survival rate 

Initial yearling survival following release 
Spawning contribution of one-ocean fish 
Fraction recovered as adults in calibration years 
CV of hbaserv 

Conditional contribution to wild escapement 
Terminal harvest rate 


Additional adult supplementation to wild 

Hatchery collection efficiency for natural fish 

Collection efficiency for hatchery fish 

Minimum natural escapement before hatchery take 

Target take of natural fish by hatchery 

Total adults needed by hatchery 

Minimum acceptable level of egg capacity 

Hatchery egg capacity 

Fraction of fry released as subyearlings 

Level of interaction between natural and released 
hatchery subyearlings (0-100) 


Fraction of -ocean fish harvested in ocean (1</<4) 
Fraction of -ocean fish harvested in mainstem river 
Fraction of -ocean fish reaching next age class 
Fraction of -ocean fish escaping to subbasin 


Mean juvenile passage survival rate 
CV of juvenile passage 


Number of games per simulation 
Number of years per game 
Selects juvenile production function 
1 = logistic response 
2 = Beverton-Holt 
3 = Ricker 
4 = empirical model (other parameters 
apply) 


Natural 
Production 
of Juveniles: 
Version 1 


natural fry can be incorporated. Alternatively, Version 2 describes the 
steps between spawning and smolt production less explicitly, relying on 
historical information on adult production to generate smolt numbers. 
Each version is described in more detail below. 

The following description applies to the base model. Specific changes in 
model structure have been made in various versions to accommodate spe- 
cial circumstances or life history traits peculiar to individual stocks, or to 
simulate specific hatchery scenarios. Such changes are easily made. 


The number of female spawners (FSPAWN) determines the number of 
naturally produced eggs. The number of female spawners is assumed to be 
binomially distributed with N = the number of spawning adults, and p = 
ffem, the probability of being female. The number of eggs is then simulated 
as a random normal deviate with mean = FSPAWN ¢egfem and standard 
deviation = FSPAWN ¢stdegg, rounded to the nearest positive integer. 
Since the number of eggs produced is the sum of reasonably large numbers 
of independent random variates (individual spawnings), the normal ap- 
proximation seems justified. 

The SLCM’s only density-dependent mechanisms are in the egg-to- 
presmolt transition. Density-dependent mechanisms operating within the 
first year or two of life are assumed to play an important role in the dy- 
namics of salmonid populations. The expected probability of survival from 
the egg stage to just before leaving the subbasin (PRESMOLT) is calcu- 
lated as a decreasing function of initial egg numbers and hatchery fry re- 
leases (details below). 

Numerous studies of the relationship between population density and 
survival indicate that the relationships—when detectable at all—are 
rather imprecise (Peterman 1987; Rothschild 1986). In addition, factors 
other than population density can contribute to large variations in early 
freshwater survival (Alderdice and others 1977; Kocik and Taylor 1987; 
Randall and others 1987). For this reason, the SLCM uses a binomial-beta 
distribution to generate the number of presmolts produced each year. The 
binomial-beta distribution increases the variance in the number of pre- 
smolts produced from a given number of eggs beyond that introduced by 
the binomial process alone. The coefficient of variation in egg-to-presmolt 
survival (cvegsv) is specified by the user as part of the input parameters. 
The model uses the expected probability of survival calculated from the 
density-dependent relationship and cvegsv to calculate the parameters of 
the beta distribution used in this transition. 

The user can choose among three options for the density-dependent rela- 
tionship used to calculate the expected egg-to-presmolt survival. The 
choices include the Beverton-Holt relationship (Beverton and Holt 1957), 


Iph 
K(p) = sl sla asl (4) 
[1+ (alphaEggs)] 
beta 
the Ricker relationship (Ricker 1954), 
E(p) = alpha[e'?#2E6¢s)] (5) 


or the logistic response function, 


(p) Wowitacl se aie ane, (6) 


1+ e(alpha+betaEGGS) 


The advantage of the logistic function is that it allows greater flexibility 
in defining the shape of the relationship between egg density and survival 
(fig. 3). With the Beverton-Holt and Ricker models, the greatest decline in 
survival always occurs at the lowest densities. With the logistic function, 
we can establish threshold levels where the density-dependent effects are 
most acute and survival decreases most rapidly. Note that the logistic 


model produces a dome-shaped egg-to-presmolt relationship similar to the 
Ricker model. 


(A) 


Ricker function 


Beverton-Holt 
function 


Survival 


Eggs. 


(B) 


Survival 


Figure 3—General shape of the 
relationship between initial egg 
density and egg-to-smolt survival 
generated using the Beverton-Holt 
and Ricker models (A) and various 
parameter combinations with the 

Eggs logistic response function (B). 


Alternative Model 
of Juvenile 
Production 
(Version 2) 


The SLCM includes the option of integrating hatchery-released fry or 
subyearlings into the simulated natural populations. In the base model, 
the number of fry released by the hatchery (FRYREL) is first adjusted for 
post release mortality using the survival parameter, hfry_sv. The surviv- 
ing fry interact with the natural fry according to the parameter, fryint. 
This parameter expresses interaction on an arbitrary scale that ranges 
from 0 to 100. If fryint = 0, then hatchery fry have no effect on survival of 
natural fry; if fryint = 100, hatchery fry affect the survival of natural fry to 
the same extent as additional natural fry. Expected survival of naturally 
spawned fish and hatchery-released fish is calculated in a multistep pro- 
cess outlined in appendix A. Admittedly, this process is arbitrary and 
rather convoluted. The method has two desired properties: (1) survival of 
naturally spawned fish decreases when hatchery-released fish are present 
and (2) the impact of hatchery releases on naturally produced fish is di- 
rectly proportional to the degree of competition. 

The base SLCM accommodates up to four year-classes of migrants. Sub- 
yearling migrants (SMOLTO) exit the tributary during the same year they 
are spawned; yearling and later migrants (SMOLT1 — SMOLTS) exit in 
succeeding years. Each additional year that the presmolts remain in the 
tributary they are subject to a mortality factor reflected in inbsmsv. 


In many instances, we do not have the stock-recruitment information 
needed to unambiguously choose one form of the juvenile production func- 
tion over other alternatives. Available data may be scarce, providing no 
hint of the underlying relationship, or the data may be so noisy that no 
single function provides a superior fit. Where few or no data are available, 
little can be done to inspire confidence in the quantitative output of any 
model. Information from similar stocks might help, but in the end the 
chosen relationship must be viewed as a hypothesis entailing considerable 
uncertainty. 

In instances where data are available but the underlying relationship is 
obscured, an empirical approach based on conditional probabilities may be 
used. This offers a way to capture the variation in recruitment, and to 
some extent, the density-dependent relationships that might exist. Data 
on spawning stock size and the estimated number of smolts produced from 
each spawning are compiled. Lognormal distributions are independently 
fit to the spawning and smolt production data; these are known as the 
marginal distributions. Smolts produced is plotted against spawning es- 
capement in a two-dimensional graph. Reference lines are drawn to sepa- 
rate the marginal distributions into thirds, thus dividing the plot into nine 
cells (fig. 4). Dividing the number of points in a given cell by the sum of 
the three cells in its column provides an estimate of the conditional prob- 
ability that the number of smolts produced in a given year will fall within 
that region of the graph, given that the corresponding spawning escape- 
ment is within the indicated range. These probabilities, along with the pa- 
rameters defining the marginal distributions, must be specified in the 
SLCM parameter file if this option is used. 

The SLCM uses these probabilities in a three-step process. First, it 
takes the number of spawners and determines which third of the spawning 
distribution is represented. Then, using the conditional probabilities, it 
randomly selects from a multinomial distribution the region of the smolt 
distribution that will be used to generate the number of smolts produced. 
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mean 
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Figure 4—Natural logs of stock-recruitment data plotted in two 
dimensions. Dashed lines separate marginal distributions of stock 
and recruits into equal thirds, allowing computation of conditional 
probabilities. 


Finally, it randomly selects a smolt production level from within the se- 
lected region of the marginal, lognormal distribution of smolts. 

The SLCM does not allow for density-dependent interactions between 
natural fry and hatchery-released fry when Version 2 is used. 


The number of eggs produced in the hatchery (HEGGS) is a random vari- 
ate drawn from a normal distribution, where the mean and variance of the 
distribution are a function of the number of females spawned in the hatch- 
ery and the specified fecundity parameters. If hatchery egg production 
does not meet the minimum level defined by htarget *eggcap, then addi- 
tional eggs (from an external source) are added to the hatchery stock. The 
number of fry produced is a binomial random variate with N = HEGGS, 
and p =egg_fry. The fraction indicated by fryplant is released as sub- 
yearlings; those remaining are released as yearling smolts the following 
year. The parameters, fry_smt and hsmt_sv, determine the survival from 
fry to yearling smolts within the hatchery and the initial survival of smolts 
shortly following release. 


The migration of smolts from upstream rearing areas to the estuary can 
be a critical period for many stocks. It is especially acute for those that 
make extensive migrations. For upper Columbia and Snake River popula- 
tions, the migration path may include up to nine mainstem hydroelectric 
dams and associated reservoirs. Thus, there are good reasons for having a 
separate component for the juvenile migration period (juvenile passage). 


Ocean and 
In-river Allocation 
of Adults 


Juvenile passage is expressed in the SLCM as a binomial-beta process. 
The number of smolts surviving downstream migration is a binomial ran- 
dom variate, where the survival probability (p) is drawn from a beta distri- 
bution. The expectation and variance of p can be estimated using auxiliary 
passage models that address this stage of the life cycle in more detail. In 
our analyses of Columbia Basin stocks, we used CRiSP, a passage model 
developed at the University of Washington, to estimate survival from 
tributary rearing areas to the Columbia River Estuary. (Information on 
CRiSP can be obtained from James Anderson, Center for Quantitative Sci- 
ences, University of Washington, Seattle, WA.) The parameters of the beta 
distribution were estimated from a set of survival estimates based on ran- 
domly generated flow conditions. In earlier analyses, we used a fourth- 
order autoregressive model with a random normal error term that was fit 
to CRiSP output. Results from both approaches were similar. Other pas- 
sage models can easily be substituted for CRiSP. 


The advent of the coded wire tag (CWT) has permitted large-scale mark 
and recapture experiments on Pacific salmon from the Columbia River Ba- 
sin. Of the millions of tagged fish released each year, however, only a 
small fraction of the tags are recovered. These recoveries provide a large 
bank of information on the age-specific relative distribution of recaptured 
adults among the several recapture sites, including ocean fisheries, in-river 
fisheries, and hatchery racks or other in-basin recovery sites. Unfortu- 
nately, the tag recovery data will never be sufficient to establish the abso- 
lute, age-specific probabilities of dying in the ocean, returning to the river, 
being harvested in the ocean or river, or dying in the river—parameters 
often included in management models. These population parameters can 
never be known precisely because there are more unknown parameters 
than recovery sites. Fish that die of natural causes are never seen. Thus, 
one cannot be sure if they died early or late in life. The best that CWT 
data can suggest is the magnitude of the total mortality loss—that propor- 
tion of the population that is never recovered or otherwise accounted for. 

Rather than modeling age-specific mortality rates and maturity sched- 
ules, the SLCM is set up to correspond to the empirical model suggested by 
CWT recovery data. The total number of adults recovered is simulated us- 
ing a binomial-beta distribution. Again, the beta distribution is used to in- 
troduce added stochasticity. A multinomial process is used to allocate the 
total number of fish recovered by age class among ocean harvest, in-river 
harvest, and return to the subbasin. Appendix B shows how the model cal- 
culates adult recovery parameters from CWT data. 

Other models that simulate ocean and in-river harvest, such as a model 
developed by the Oregon Department of Fish and Wildlife (1991), can be 
used to estimate ocean and in-river recovery probabilities. In such cases, 
the simulated harvest and escapement data are used instead of CWT data. 
Harvest patterns that change with time are accommodated by using a spe- 
cial subroutine that calculates cohort-specific recovery probabilities using 
simulation results from the harvest model. Since the total fraction recov- 
ered depends on harvest and mortality rates, adtrecv is adjusted for each 
cohort to reflect changing conditions. 


Subbasin Allocation of adults returning to the subbasin follows a multinomial pro- 
Allocation cess, with the twist that hatchery needs and natural escapement targets 
have priority over harvest. Hatchery and natural adults are allocated 
among hatchery spawning, natural spawning, and terminal harvest accord- 
ing to a suite of parameters that direct the allocation process shown in fig- 
ure 5. The number of fish available to the hatchery is determined by the 


Hatchery 
Returning Adults 
(hsubesp) 


Natural 
Returning Adults 


(subesp) 


<——hcolrate 


Hatchery Trap 


minesp 
maxthar Ip ce 


subharv 


prspsv ———— <- ———— hprspsv 


hsubharv 


Natural Spawning 


Hatchery Spawning 


A 
I 


prspsv 


Figure 5—Conceptual flow diagram of hatchery take and natural 
escapement of adults. The abbreviations are explained in table 1. 
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Operational Notes 


hatchery trap efficiency, which can be different for natural and hatchery 
fish. Of the natural fish that are available, the hatchery removes the 
smaller of hatrem (the hatchery target) or those available in excess of 
minrun (the minimum acceptable wild escapement). The remainder of 
hatchery needs are met using returning hatchery adults, if available. 
Adults not needed by the hatchery and those escaping the trap are subject 
to a terminal harvest, with independent rates for hatchery and natural 
fish. Hatchery adults which escape harvest stray into natural spawning 
areas in proportion to the parameter, hstray. (This parameter may also in- 
corporate the differential spawning success between natural and hatchery 
fish in the wild.) The total number of adults spawning naturally includes 
both natural and hatchery adults, and an additional number of adults, 
adtsup, which originate from outside the system (simulating adult supple- 
mentation). All adults are subject to prespawning mortality, either within 
the hatchery or natural areas. Jacks (one-ocean fish) can be excluded from 
the spawning population by using the parameters, jackspn and hjackspn. 


The SLCM and ancillary programs are written in the SAS programming 
language. While the model has been developed on a SUN workstation, it 
can be transferred easily to any system equipped to run the SAS language. 
Copies of files needed to run the model and produce summary analyses are 
available from the authors. The modular code can be modified to suit the 
particular needs of various user groups. 

Once the SLCM files are installed, users run the model by setting up 
three parameter files, a file containing initial numbers, and a control file. 
Examples are shown in appendix C. These files can be modified by using 
a convenient text editor. Only the control file needs to be submitted as an 
SAS program. Ifthe control file is set up properly, it will make all neces- 
sary calls to other files and programs. As output, the model creates an 
SAS data set containing all variables requested in the control file. Model 
results can be analyzed using ancillary SAS programs supplied with the 
model. Users may write additional ancillary programs to meet their spe- 
cific needs. 

The control file includes macros that allow users to perform multiple 
runs of the model using stepped parameters (each run would consist of 
multiple games), and to change parameters for specific years within games. 
This is useful in performing sensitivity analyses and in simulating environ- 
mental changes over time. The stochastic properties of the model can be 
changed through this mechanism by introducing additional uncertainty or 
changing the distributional properties of the parameters. 

The control file also contains an SAS macro variable, “HATCHERY. 
Setting “HATCHERY to “off’ shuts down all hatchery production. The 
simulated population then will contain only naturally produced fish. 


ILLUSTRATIVE MODEL RESULTS 


We demonstrate the SLCM using parameter sets listed in appendix C. 
These parameters were chosen to reflect a stock that is not limited by habi- 
tat, but by a low smolt-to-adult survival rate. Thus, the density-dependent 
parameter (beta) is set at a very low level. In this example, subbasin es- 
capement (natural and hatchery) is the principal output examined. Our 
base case is a contrived fall chinook example that is not meant to mimic 


any particular stock. The base case includes an assumption of little inter- 
action between hatchery and natural fish, small amounts of fry releases 
(the bulk of hatchery releases are smolts), and no adult supplementation. 
Almost all surplus hatchery fish are harvested. The logistic response func- 
tion is used to describe juvenile production. 


Single Game The simulation of a single population over a 500-year period is shown 
in figure 6. The population’s fluctuation over time is typical of the fluctua- 
tions of state variables within the SLCM. Though we have not rigorously 
compared the model output with empirical time series, the results appear 
consistent with general observations of populations in nature. 

A second way of viewing simulation results is to look at the relationship 
between the number of adults spawning naturally in a given year and their 
offspring recovered in future years as adults. In figure 7, SPAWNERS and 
associated RECRUITS for the 500-year period are plotted along with deter- 
ministic representations of the underlying spawner-recruit relationship. 
The curved line represents the density-dependent relationship between the 
number of spawning adults and the expected number of recruits. The 
straight line is the inverse of the relationship between a given number of 
recruits and the number of spawners they would produce. In a determinis- 
tic model, the population would move quickly to the equilibrium point de- 
termined by the intersection of the two lines. 

In contrast, the SLCM suggests no single equilibrium point, but rather 
a broad clustering of points located near the deterministic equilibrium 
points. The first moment of the mass of point scattered within the spawner- 
recruit plane is actually located below and to the left of the intersection 
point. This suggests that while stochastic processes in the SLCM occasion- 
ally lead to production levels much higher than average, average produc- 
tion is less than predicted by deterministic models. 


Year 


Figure 6—Subbasin escapement for a single simulation over time. Solid 
line indicates natural fish, dashed line indicates hatchery fish. 
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Figure 7—Simulated data showing the number of adult recruits 
produced by a given number of spawning adult salmon. The curved 
line represents the density-dependent relationship between the 
number of spawning adults and the expected number of recruits. 
The straight line is the inverse of the relationship between a given 
number of recruits and the number of spawners they would produce. 


Since the stochastic nature of the SLCM produces a different outcome 
for every game or simulation, the outcome from a single game is of little 
value. It is more informative to run multiple games and examine the out- 
put collectively. In this manner, one can evaluate possible future states 
probabilistically. For example, we ran 500 games covering 50 years each 
using the SLCM. When combined, the output from these games creates a 
probability distribution for each year, such as the distribution of spawners 
in year 50 (fig. 8). The distribution shown in figure 8 is characteristic of 
the positively skewed distributions of state variables common to SLCM 
results. 

Portraying the distribution of a state variable through time is problem- 
atic. In figure 9, contour lines show the shape of the distribution of spawn- 
ers through the simulation period. The contour lines depict the minimum, 
10th percentile, median, mean, 90th percentile, and maximum from the 
500 games. For this example, the bulk of the distribution is fairly stable, 
with the most variation occurring in the upper tail of the distribution. 
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Figure 8—Distribution of the number of natural spawners in year 50 
from 500 simulations. 
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Figure 9—Distribution of natural spawners over time from 500 simulations. 
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Paired 
Comparisons 


Models are often used to examine the relative effectiveness of alternative 
management proposals for increasing fish production. The effects of man- 
agement changes could be lost or confounded with random fluctuations in- 
troduced by stochastic models such as the SLCM. To facilitate meaningful 
comparisons of alternative management proposals, the SLCM is designed 
to allow the user to control the random number streams used within the 
model. This allows users to replicate serial patterns within the SLCM re- 
flecting changes in environmental conditions over time. Thus, a user could 
compare two or more management alternatives over the same temporal 
patterns in survival. For example, figure 10 shows two time series for a 
matched pair of games where downstream migrant survival has been in- 
creased by 40 percent in the second game. While there is independent 
variation in each game, a significant portion of the variation is shared. 
Multiple paired comparisons can lead to distributions of comparison statis- 
tics in much the same way that multiple games lead to distributions of 
state variables. 


SENSITIVITY ANALYSIS 


The structure of the SLCM allows explicit inclusion of parameter uncer- 
tainty in the model. By increasing the variance of survival processes, users 
can include some of the existing uncertainty about survival or production 
parameters. But increasing the variance only addresses questions of preci- 
sion; it does not directly address problems of accuracy that arise if the pa- 
rameter estimates are biased. A rigorous sensitivity analysis is needed to 
assess the ramifications of potential bias in parameter estimates. 

The SLCM’s relatively straightforward structure allows sensitivity 
analyses to be performed easily. The model can be reduced to five major 
components: density-independent juvenile production and survival, 
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Figure 10—Paired comparison of model results using two levels 
of passage survival: 32 percent (lower line) and 44.8 percent 
(upper line). 


Results of the 
Sensitivity 
Analysis 


density-dependent juvenile survival, downstream passage, adult recruit- 
ment, and harvest. The expected system response in each component can 
be reduced to a linear combination of certain model parameters (or a single 
parameter in some cases). Thus, users can select a small subset of param- 
eters to serve as surrogates for all parameters when conducting a sensitiv- 
ity analysis. 

In a preliminary sensitivity analysis, we examined the sensitivity of the 
natural population to parameter changes by systematically varying logitscl 
and adtrecv to reflect different levels of early juvenile and postsmolt sur- 
vival, varying mpass to reflect changes in downstream migrant survival 
(passage survival), and varying the fraction of recruits within a given age 
class that escapes to the subbasin (as opposed to being harvested). These 
four factors were varied in combinations to create two-dimensional matri- 
ces of parameter sets (either 11 x 11 or 10 x 10, depending on the vari- 
ables). Within each matrix, logitscl, mpass, and adtrecv varied by a factor 
of three, while the proportion of nonharvested fish varied by a factor of ten 
(0.1 to 1). Five hundred games were simulated for each parameter set. 
Number of spawners in years 15 and 100 were compared across all param- 
eter combinations. The number of spawners in year 15 proved to be a good 
indicator of short-term stock dynamics, while the 100-year results were 
better indicators of long-term effects. 


Some comments on the sensitivity of the SLCM to changes in parameters 
are appropriate, though we refrain from a detailed discussion. First, the 
distribution in number of spawners is very sensitive to changes in param- 
eter values of the magnitude that we examined. Most sensitive is the up- 
per tail of the distribution, which changes rapidly following parameter 
changes. 

Second, short-term model results are more sensitive to changes in adult 
survival than to changes in juvenile survival. The mean number of spawn- 
ers in year 15 was greater following increases in adult survival than after 
equivalent increases in juvenile survival. This is understandable because 
of the time lag required before changes in juvenile survival manifest them- 
selves in changes in spawning escapement. Given a longer time period 
(such as 100 years), equivalent changes in juvenile survival and adult sur- 
vival produced equivalent effects. However, the management implications 
of the apparent short-term effectiveness of increasing adult survival should 
not be ignored. 

Among the variables that we examined, the model first appeared to be 
most sensitive to changes in passage survival. On further inspection, the 
apparent sensitivity to passage survival could be explained by the rela- 
tively low variance used for passage survival compared to the variance 
used for juvenile and postsmolt survival. The SLCM, like other stochastic 
models of this type, is most sensitive to the survival factor with the lowest 
variance. The interplay between survival and variance in survival can be 
demonstrated by plotting the changes in the mean number of spawners as 
the coefficient of variation (CV) in egg-to-adult survival increases (fig. 11). 
While the decline in the mean number of spawners is significant, the in- 
crease in the CV of state variables such as the number of spawners is more 
dramatic. We observed that the relationship between the CV of spawners 
and the CV of survival had a slope that is very close to (but slightly greater 
than) one. 
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Figure 11—Effect of the coefficient of variation (CV) of egg-to- 
adult survival on the mean number of spawners in year 30 and 
the CV of spawners. 


MODEL CALIBRATION 


Proper calibration of a model such as the SLCM is always problematic. 
Rarely, if ever, is there sufficient information to provide independent esti- 
mates of all model parameters. For many stocks, the most readily avail- 
able information likely will be that collected from hatchery fish, such as 
sex ratios and fecundity information. Given the large number of hatcher- 
ies within the Pacific Northwest, reliable information on most stocks or on 
stocks similar to those in question should be available. Users should be 
able to obtain fairly accurate estimates of the average potential egg produc- 
tion from a given number of spawners and estimate a variance. Howell 
and others (1985) and the recently completed subbasin production plans 
(Columbia Basin Fish and Wildlife Authority 1990) are useful sources of 
this type of information. 

The second most readily available source of information is CWT data. 
The Pacific Marine Fisheries Commission maintains a comprehensive data 
base on tagged fish releases and recoveries for Columbia River and other 
Northwestern stocks. Further information on using CWT data to estimate 
harvest and escapement parameters is provided in appendix B. Two as- 
pects of CWT data should be remembered: (1) most CWT data are from 
hatchery fish, and (2) the ratio of tags returned to the total number of 
tagged fish released generally is extremely small. Because most data are 
from hatchery fish, caution is called for when extrapolating data to natu- 
rally produced fish. Because the number of returns is so small, acute bi- 
ases in parameter estimates might arise because of uneven reporting rates 
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A General Process 
for Calibration 


from the different fisheries or recovery sites. In general, CWT data do not 
provide reliable estimates of the total recovery fraction (adtrecv) because 
of low, uncertain sampling rates. They likely do provide, however, the best 
available indication of the allocation of recovered fish among fisheries and 
subbasin returns. 

Reliable estimates of survival parameters are very difficult to obtain. In 
an ideal situation, modelers would have independent data on spawner den- 
sities and smolt production, juvenile passage survival, and survival from 
the estuary through adulthood. Such data do not exist. For most stocks, 
the best available data will be in the form of a continuous record or time 
series of counts of individuals at a given life stage, such as estimates of the 
number of adults returning to the confluence of a major tributary. The 
trick to calibrating the model is to use these data in combination with all 
other sources of information to derive a set of parameter estimates that 
provides the best fit to the data. 


Calibration depends in part on the version of juvenile production that is 
used. When using the more conventional model of juvenile survival (Ver- 
sion 1) is used, we found the following process to be successful. First, we 
obtained estimates of as many parameters as possible from the literature, 
opinions of experts, other models, and CWT data. From this information, 
we estimated all but the juvenile survival parameters and total adult re- 
covery fraction. 

Second, we estimated the general shape and magnitude of the egg-to- 
presmolt survival function. Direct estimates of smolt production resulting 
from different spawning levels generally are not available. Even if stock- 
recruit data are available, they usually are too scattered to allow precise 
estimation of a survival function. Ultimately, considerable professional 
judgment is needed to construct an approximate juvenile production func- 
tion. The term “production function” refers to S¢f(S), where S is the num- 
ber of eggs, and f(S) is the egg-to-presmolt survival function. Since the pa- 
rameters of the production function and survival function are the same, the 
terms are interchangeable at times. 

The principal factors to consider in constructing a juvenile production 
function are: (1) Does the function have the proper shape? (2) If the pro- 
duction function has a maximum, does it correspond to the level of spawn- 
ing adults that would be expected to produce the maximum number of 
presmolts? (3) Do the relative differences in egg-to-presmolt survival at 
low density and at high density seem reasonable? 

It is critical that the shape of the production function for juveniles be cor- 
rect, since this function directly determines the compensatory capacity of 
the population. Accurate estimates of the absolute number of presmolts 
produced are less important because our principal interest is in the number 
of adults produced. In the final step (described below), we scale the adtrecv 
parameter so that the combined parameter set provides the best fit to the 
available data, and compensates for some of the error we might have intro- 
duced through estimates of other parameters. 

The final step in calibration is to estimate adtrecv, cvegsv, and cvadtrv. 
For this purpose, a second ancillary model has been constructed. The cali- 
bration model uses the parameters of the SLCM and deterministic analogs 
of its relationships to predict elements of a historic time series using previ- 
ous observations in the time series. The parameter adtrecv is introduced 
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Calibration Using 
the Alternative 
Model of Juvenile 
Production 


as a vector of annual coefficients that is adjusted systematically to mini- 
mize the sum of squared errors between observed and predicted observa- 
tions. Users can set constraints on the annual coefficients so that they fall 
within a given range. The mean of the fitted vector is then used as an esti- 
mate of adtrecuv. 

The calibration model also calculates the mean and variance of the prod- 
uct of egg-to-presmolt, downstream migrant, and smolt-to-adult survival 
for each cohort. The variance of this product is partitioned among the 
three survival components algebraically. Since there is only one variance 
equation but three unknowns, the user must decide the relative propor- 
tions of the variance to be allocated to each life stage, or alternatively must 
set one or two of the component variances to a constant and solve for the 
remainder. 


The calibration process for the alternative model of juvenile production 
(Version 2) is similar in concept to that described above. Smolt-to-adult 
survival again is used as a scaling factor to provide the best fit to the data. 
Paulsen and others (1991) provide an example using this version of the 
model. 


CONCLUSIONS 


The SLCM has proven to be a powerful and flexible tool for examining 
the impact of changes in population parameters on the structure and size 
of salmonid populations. While some aspects of the SLCM’s behavior have 
been difficult to understand at first glance, the model’s overall behavior 
conforms well to theory and experience with deterministic models. Some 
of the more provocative results are those that deviate from an analogous 
deterministic model, such as the model’s sensitivity to parameter variance. 
As our experience with the model grows, we expect new questions to arise 
concerning the interplay of stochastic factors and population dynamics, 
opening new areas of research. Such questions, and their answers, could 
have significant management implications. 


REFERENCES 


Alderdice, D. F.; Bams, R. A.; Velsen, F. P. J. 1977. Factors affecting 
deposition, development, and survival of salmonid eggs and alevins: a 
bibliography, 1965-1975. Tech. Rep. 743. Nanaimo, BC: Canada Fisheries 
and Marine Service, Pacific Biological Station. 276 p. 

Beverton, R. J. H.; Holt, S. J. 1957. On the dynamics of exploited fish 
populations. Fishery Investigations, Ser. 2, Vol. 19. London, UK: 
Ministry of Agriculture, Fish and Food. 533 p. 

Boswell, M. T.; Ord, J. K; Patil, G. P. 1979. Chance mechanisms underlying 
statistical distributions. In: Ord, J. K.; Patil, G. P.; Tallie, C., eds. 
Statistical distributions in ecological work. Fairland, MD: International 
Co-operative Publishing House: 3-156. 

Columbia Basin Fish and Wildlife Authority. 1990. Columbia Basin system 
planning salmon and steelhead production plans. Portland, OR: North- 
west Power Planning Council: 31 volumes. 

Gilpin, M. E.; Soulé, M. E. 1986. Minimum viable populations: processes of 
species extinction. In: Soulé, M. E., ed. Conservation biology. Sunderland, 
MA: Sinauer Associates: 19-34. 


22 


Goodyear, C. P. 1977. Assessing the impact of power plant mortality on 
the compensatory reserve of fish populations. In: Van Winkle, W., ed. 
Proceedings of the conference on assessing the effects of power-plant- 
induced mortality on fish populations. New York, NY: Pergamon Press: 
186-195. 

Howell, P.; Jones, K.; Scarnecchia, D.; [and others]. 1985. Stock assessment 
of Columbia River anadromous salmonids: Vols. I and II. DOE/BP-12737- 
1,2. Portland, OR: U. S. Department of Energy, Bonneville Power 
Administration, Department of Fish and Wildlife. 1032 p. 

Jensen, A. L.; Hamilton, T. A. 1982. Application of a conventional fishery 
model for assessment of entrainment and impingement impact. 
Environmental Biology of Fishes. 7: 181-185. 

Kocik, J. F.; Taylor, W. W. 1987. Effect of fall and winter instream flow on 
year-class strength of Pacific salmon evolutionarily adapted to early fry 
outmigration: a Great Lakes perspective. In: Dadswell, M. J.; Klauda, 

R. J.; Moffitt, C. M.; [and others], eds. Common strategies of anadramous 
and catadromous fishes. American Fisheries Society Symposium. 
1: 430-440. 

MacCall, A. D.; Parker, K. R.; Leithiser, R.; Jessee, B. 1983. Power plant 
impact assessment: a simple fishery production model approach. Fishery 
Bulletin. 81: 613-619. 

McConnaha, W. E. 1992. System Planning Model documentation update 
version 5.10. Portland, OR: Northwest Power Planning Council. 28 p. 

Oregon Department of Fish and Wildlife. 1991. Preliminary evaluation 
of the Northwest Power Planning Council’s proposed amendments for 
Snake River fall chinook. Written comment submitted to Northwest 
Power Planning Council, October 16. Portland, OR: Oregon Department 
of Fish and Wildlife. 

Paulsen, C. M.; Lee, D. C.; Hyman, J. B.; Wernstedt, K. 1991. Mid- 
Columbia propagation and passage cost-effectiveness analysis. Report to 
the Bonneville Power Adminstration, September 17. Washington, DC: 
Resources for the Future. 65 p. 

Peterman, R. M. 1982. Model of salmon age structure and its use in 
preseason forecasting and studies of marine survival. Canadian Journal 
of Fisheries and Aquatic Sciences. 39: 1444-1452. 

Peterman, R. M. 1987. Review of the components of recruitment of Pacific 
salmon. In: Dadswell, M. J.; Klauda, R. J.; Moffitt, C. M.; [and others], 
eds. Common strategies of anadramous and catadromous fishes. 
American Fisheries Society Symposium. 1: 417-429. 

Pielou, E. C. 1977. Mathematical ecology. New York: John Wiley. 385 p. 

Press, S. J. 1989. Bayesian statistics. New York: John Wiley. 237 p. 

Randall, R. G.; Healey, M. C.; Dempson, J. B. 1987. Variability in length 
of freshwater residence of salmon, trout, and char. In: Dadswell, M. J.; 
Klauda, R. J.; Moffitt, C. M.; [and others], eds. Common strategies of 
anadramous and catadromous fishes. American Fisheries Society 
Symposium. 1: 27-41. 

Ricker, W. E. 1954. Stock and recruitment. Journal of the Fisheries 
Research Board of Canada. 11: 559-623. 

Rothschild, B. J. 1986. Dynamics of marine fish populations. Cambridge, 
MA: Harvard Press. 277 p. 


23 


Taylor, M. W. 1981. A generalized inland fishery simulator for 
management biologists. North American Journal of Fisheries 
Management. 1: 60-72. 

Walters, C. J. 1981. Optimum escapements in the face of alternative 
recruitment hypotheses. Canadian Journal of Fisheries and Aquatic 
Sciences. 38: 678-689. 

Walters, C. J. 1986. Adaptive management of renewable resources. New 
York: Macmillan. 374 p. 


24 


I. Definitions 


II. Calculations 


APPENDIX A—A PROCESS FOR ADJUSTING EGG-TO- 
PRESMOLT SURVIVAL WHEN HATCHERY FISH ARE 
INTRODUCED 


pl = probability of survival of naturally produced eggs with no additional 
hatchery fish 


p2 = pl adjusted for presence of hatchery fish 
p3 = p2 adjusted for time period 

N(¢) = natural fry alive at time t 

H(t) = hatchery fry alive at time t¢ 

u1 = instantaneous mortality rate 

u2 = adjusted mortality rate 

E(¢) = denotes expectation 


Binomial{N,p} = function call that generates a random number from 
a binomial distribution. 


Beta{p,cv} = function that generates a random number from a beta 
distribution 


f{x} = function call to Genet dependent relationship that defines p 
as a function of x 


All other definitions are as aidefned within the text. 


Calculations take place in sequence 


. Ll = -log(p1)/100 

. N(reltime) = Binomial{EGGS,exp(-11*reltime) 

. H(reltime) = Binomial {hatfry,hintsuv 

. E(p2) = f{(EGGS+H(reltime)} 

. "2 = —log(p2)/100 

. E(p3) = exp(—u2*(100-reltime)) 

. pd = Beta{E(p3),cvegsv} 

. PRESMOLT = Binomial{N(reltime)+H(reltime),p3} 


ONonhwaoNdre 
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APPENDIX B—AN EXAMPLE OF HOW TO ESTIMATE 
OCEAN AND IN-RIVER PARAMETERS FROM CWT DATA 


Snake River Fall Number released: release = 1,467,443 


Chinook CWT 
Data Recoveries: 
Ocean In-river Hatchery 
harvest harvest returns 
ol = 2,105 rl = 1,390 hl = 1,078 
o2 = 12,608 r2 = 4,716 h2 = 188 
03 =2, 137. r3 = 1,938 h3 = 142 
04 = 92 r4 = 121 h4 =4 
Assumption Survival from release site to below Bonneville dam: surv = 0.3 
Calculations total =01+02+03+044+r1+r2+r3+r4+hl+h2+h3+h4 
suml =total—(ol+rl1+h1l) 
sum2 =suml—(02+r2+h2) 
sum3 =sum2-—-—(03+r3 +h3) 
adtrecv = total / (release ® surv) adtrecv = 0.06 
ocnlhav = o1/ total ocnlhav = 0.079 
ocnisuv = sum1 / total ocnisuv = 0.828 
ocnirvh = rl / total ocnirvh = 0.052 
ocnlsbe = hl/ total ocnlsbe = 0.041 
ocn2hav = 02/suml ocn2hav = 0.575 
ocn2suv = sum2/suml1 ocn2suv = 0.202 
ocn2rvh = r2/suml ocn2rvh = 0.215 
ocn2sbe = h2/sum1 ocn2sbe = 0.009 
ocnshav = 03 / sum2 ocnshav = 0.482 
ocn3suv = sum3 /sum2 ocn3suv = 0.049 
ocn8rvh = r3/sum2 ocn38rvh = 0.437 
ocn3sbe = h3 / sum2 ocn3sbe = 0.032 
ocn4hav = 04/ sum3 ocn4hav = 0.424 
ocn4rvh = r4/sum3 ocn4rvh = 0.558 
ocn4sbe = h4/sum3 ocn4sbe = 0.018 
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C.1 Control File 


APPENDIX C—EXAMPLE CONTROL FILE AND PARAMETER 
INPUT FILES 


KKK KKKKEKKEKKE KKK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKKKKKKK KKK KKKKEKEK 


* CONTROL FILE FOR RUNNING THE STOCHASTIC LIFE CYCLE MODEL * 
KEKKKKKKKKKKKKKKKKK KKK KK KKK KKK KK KKK KK KKK KK KKK KKK KKK KKK KKK KKK KEK KK KK KK KKK KKKKEKEK 
* CHANGES IN THIS FILE SHOULD BE MADE ONLY TO LOWER CASE IDENTIFIERS. * 
* DO NOT DELETE ANY SEMICOLONS OR OTHER SPECIAL CHARACTERS! ! ! * 


KKK KKK KKK KK KK KEK KK KK KEK KKK KKK KKK KK KKK KKK KKK KK KEK KKK KKK KEKE KKK KKK KKKKKK KKK KKK KKKKK, 
’ 


%INCLUDE ‘macro6_1.sas’ /source2; 
KRKKK KK KEKE KKK KKK KEK KKK KEKE KKK KKK KKK KKK KKK KKKKKEEK 


* PART 1. CONTROL MACRO AND KEY FILES * 

RRA KAKA KKK RRK EK AK KEKE KER KKK KKK ERKEKE RK ERE EK» 

LIBNAME SASLIB ‘slcmout’; /* DIRECTORY FOR SAS OUTPUT FILES */ 
FILENAME J_PARMS ‘test61j.par’; /* FILE WITH PRODUCTION PARAMETERS */ 
FILENAME P_PARMS ‘testp.par’; /* PASSAGE PARAMETERS FOR WILD SMOLTS*/ 
FILENAME HP_PARM ‘testp.par’; /* PASSAGE PARAMETERS FOR HATCHERY SMOLTS*/ 
FILENAME A_PARMS ‘testa.par’; /* FILE WITH ADULT PARAMETERS (WILD) */ 
FILENAME HA_PARM ‘testa.par’; /* FILE WITH ADULT PARAMETERS (HATCHERY) */ 
FILENAME INITIAL ‘initnum.exp’; /* FILE WITH INITIAL NO. (WILD HAT.) */ 
SLET RESULTS = explél ; /* SAS DATASET TO CONTAIN MODEL OUTPUT */ 

SLET JUVPARM = facbase ; /* HANDLE FOR ACCESSING PRODUCTION PARAMETERS */ 
SLET PASS = $STR(‘testl’) ; /* HANDLE FOR PASSAGE PARAMETERS */ 


%LET ADULT = %STR(‘testl’) ; 
/* HANDLE FOR ADULT RECOVERIES */ 
%LET HATCHERY = %STR(‘on’) ; /* TURNS HATCHERY PRODUCTION ON OR OFF */ 


SMACRO CONTROL; 


NGAMES = 500; /* NUMBER OF GAMES */ 
YEARS =950;; /* NUMBER OF YEARS PER GAME */ 
SMTTOG = 1; /* SPECIFY SMOLT PRODUCTION OPTION */ 


SMEND CONTROL; 


KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KK KKK KKKKEKKKKKKKKKKKKKKK 


* PART 2. IDENTIFY STEPPED VARIABLE AND LEVELS * 
KKK KR RK KR RK RK RK KR KEK KK KKK KKK RK KK KKK KK KKK KK KKK KK EK KKK KEKE - 
%LET STEPVAR = noparm; /* DEFINE STEPPED VARIABLE */ 
SMACRO STEPDEF; 
NSTEPS = 
LEVEL1 = 
LEVEL2 = 
LEVEL3 = 
LEVEL4 = 
LEVEL5 = 
LEVEL6 = 
LEVEL7 = 
LEVEL8 = 
LEVEL9Y = 
LEVEL10 = 0 ; 
$MEND STEPDEF; 


; /* DEFINE NUMBER OF STEPS, 1 OR MORE */ 
i /* DEFINE VALUES FOR EACH STEP. i / 
i /* UP TO TEN LEVELS ARE PERMISSIBLE */ 


(oy (Oy (oy (2) (ooo) oyna 
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C.2 Production 
Parameters 


KEK KKK KKK KKK KE KKK KEK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KK KKK KEK KEEEKK 


* PART 3. IDENTIFY YEAR-DEPENDENT PARAMETERS AND VALUES * 
oe FOR VARIABLES OTHER THAN PASSAGE AND OCEAN HARVEST * 
EPR LL SILI AI AILS SIE IE TITS LILI I MII Feed Ud ES tad thea Ld Pt td Utada ESL ILI RUN FES Lata ed IER EI ESE 
%SMACRO VARYEAR; /* OPTIONAL SUBROUTINE */ 
* SELECT; 
*WHEN (0 < YEAR <= 2) DO; 
*parm = value; 
*END; 
*WHEN (2 < YEAR <= 5) DO; 
*parm = value 
*END; 
* et cetera; 
*OTHERWISE parm = value; 
*END; 
SMEND VARYEAR; 


KRKKKE KKK KKK KEK KEKE KEKE KKK KK KEK KKK KEKE KEKE KEKE KEKEKKKKKEKKEKKKKKKKEK 


* PART 4. IDENTIFY OUTPUT VARIABLES OF INTEREST = 
CFE St TEI I Mad Fad FL CLS I TAAL IST IE TE ISIE FUSES CHIE ICICI TFL Nad Sad EAE AS Ud ad aed LaF EPCS 
SMACRO KEEPLIST; /* THESE VARIABLES WILL BE INCLUDED IN OUTPUT FILE */ 
KEEP STEP &STEPVAR STOCK GAME YEAR STRATID P_OPTION H_OPTION 
subesp subharv recruits spawners hsubesp hsubharv hrecruit hspawn 
espwild; 
stratid = ‘test case’; /* IDENTIFY STRATEGY */ 
P_OPTION = &PASS; 
H_OPTION = &ADULT; 
%MEND KEEPLIST; 


KEK KKK KK KKK KKK KKK KKK KKK KKK KEKE KKK KKK KEKE KKK KKK KKEKKKKKKKKEEK 


* PART 5. OUTPUT CONTROL x 


FOCI III ICI III IO ICI IK 
SMACRO YEAROUT; 
ALLYEARS= ‘Y’ 8 /* PUT Y IF OUTPUT WANTED EACH YEAR, 
/* PUT N IF OUTPUT ONLY ON LAST YEAR 
SMEND YEAROUT; 


KKK KK KKK KKK KKK KEKE KKK KEK KK KEKE KKK KK KKK KE KKK KKK KKK KKEKKKKKEKEEE 


* PART 6. CALLS TO OTHER PROGRAMS bd 


KKKKEKKKK KER REE KR KR KE KEE KER KER KEKE KEK KEE KEK K ERE KEKE EEE REKEEEK » 


%INCLUDE ‘slcm6_1.mod’ /SOURCE2; 


xxx filename: test61j.par ; 
FACBASE: /* Identifies a specific parameter set */ 


**x*x Natural Production Parameters ***; 


stock = 3); /* stock identifier */ 

ffem = 0.583; /* fraction female */ 

prspsv = 0.85; /* prespawning survival */ 

egfem = 4297; /* mean eggs per female */ 

stdegg = 233; /* std. error of egfem */ 

jackspn = 0; /* spawning contribution of 1-ocean fish */ 
logtscl = lO}; /* logit scale parameter */ 

alpha = =1.73; /* density-indp. juvenile survival parm. */ 
beta = -8.41e-08; /* density-dep. juvenile survival parm. */ 
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Af 
sai 


C.3 Passage 
Parameter File 


cvegsv = 
stayl = 
stay2 = 
stay3 = 
inbsmsv = 
basercv = 
cvadtrv 


minesp = 
maxthar 


/* CV of egg-to-presmolt survival */ 

/* cond. prob. of staying over 1st winter */ 
/* cond. prob. of staying over 2nd winter */ 
/* cond. prob. of staying over 3rd winter */ 
/* inbasin overwinter survival */ 

/* fraction recovered as adults */ 

/* CV of basercyv */ 

/* minimum escapement goal */ 

/* maximum terminal harvest */ 


*** Hatchery Stock Parameters ***; 


hffem = 
hprspsv = 
hegfem = 
hstdegg = 
egg_fry = 
hfry_sv = 
fry_smt = 
hsmt_sv = 
hjackspn 
hbasercv 
hevadtrv 
hstray = 
htermhar 


0.583; 
0.85; 
4297; 

233; 


Se) (Oe 
OO wWwrP OI © DO 


(ey (2) oyXKS) 


“ese eee 


=e 


/* fraction female */ 

/* prespawning survival */ 

/* mean eggs per female */ 

/* std. error of egfem */ 

/* egg-to-fry survival rate */ 

/* initial subyearling surv. following release */ 
/* fry-to-smolt survival rate */ 

/* initial yearling surv. folowing release */ 

/* spawning contribution of 1l-ocean fish */ 

/* fraction recovered as adults */ 

/* CV of basercv */ 

/* conditional contribution to wild escapement */ 
/* terminal harvest rate */ 


*** Hatchery Management Parameters ***; 


adtsup = 
colrate = 
hcolrate 
minrun = 
hatrem = 
hatneed = 
htarget = 
eggcap = 
fryplant 
fryint = 


RETURN; 


filename: testp.par 


0; 

OSL; 
0.8; 
500; 
20; 
200; 
ORS 
500000; 
0.05; 


10; 


/* additional adult supplementation to wild */ 

/* hatchery collection efficiency of wild fish */ 

/* collection efficiency of hatchery fish */ 

/* minimum wild subesp before hatchery take */ 

/* target take by hatchery of wild adults */ 

/* total number of adults needed by hatchery */ 

/* minimum acceptable frac. of hat. egg capacity */ 

/* hatchery egg capacity */ 

/* fraction of hat. fry released as subyearlings */ 

/* level of interaction between natural and 
supplemented hatchery subyearlings (0 - 100) */ 


JUVENILE PASSAGE TEMPLATE 


The basic format for this data set is to place an asterisk in the first 
column of the line containing the scenario identifier followed by a 
variable which indicates how many lines of data follow that are 
specific to that scenario. Data lines consist of three variables: 

msurv (mean survival), passcv (CV of passage survival), and 

nyr (the number of years for which the passage parameters apply). 


msurv passcv nyr 


* base 5 
0.24 0.2 il 
0.26 0.21 2 
0.27 OLS 2 
53) 0.3 4 
OF32 0.31 2 
x trestles 
0.24 OR2 ab 
0.27 0.23 2 
0.3 O25 3 
O32 OF33 4 
OF35 0.4 2 
C.4 Adult filename: testa.par 


Recovery File 


ADULT RECOVERY PROFILES TEMPLATE (CONTRIVED EXAMPLE) 


The basic format places an * in the first column of the row containing 
the scenario identifer, followed by the scenario id and nk, a parameter 
which tells how many years of data follow. The first year of data MUST 
be data representing the calibration years. The general format for the 
data lines is as follows: 


begin - escapement — ocean harvest river harvest 
year number el 62) e3)) 7e4) oll 202) oS O4m irae 2 iors ae 


* base 5 
1989 LEO OME LS! 201 S21 ON S54 53.355 eso 4 6melii 22a Sea A 
1990 1000 8 Sorns5 yp USS as) aun ale Saat a7 alts) aS) 
1991 28005 125) 273 edkli2 la 49M 614 AS3 ie VAS 62923)5 33.9) 6! 
1992 LS OOM Or SIN Ap aer/2 9° 32: 389° 304 30) 40/250) 218) 39 
1993 1200 Z Th Se ay aly “algal alsa) als} alg? Gar Sig aL) 


* SEES tS 
1989 1600 19) .:20il). 282) S10, 2 35)453. 355.36 46 av2u24eme 44 
1990 LOO OM aa 7 SNS Ey ale CHOY aS: AG al a7 ally al) 
1991 2800. 385 290 1195 Sy S25655 (502) Si 42) 2505360) 65 
1992 1800 S18 LTO HS) MONE sseaZerssi 3a eae os 23 mas 
1993 1200 8 90 37 5 ale dey alsyy7 al, ak 7/7 alls) AO) 
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__ Astochastic model simulating the life cycle of anadromous salmonids that mimics the — 
basic mechanisms regulating populations of Pacific salmon is presented. The SLCM 
incorporates stochastic or probabilistic processes at each step in the life cycle. The model’s 
predictions are expressed in probabilistic terms. The model is designed to use inputs from 
more detailed models for specific life stages, in combination with a minimum number of 
empirically based parameters. Users can choose among three density-dependent relation- 
ships to describe egg-to-smolt survival and among a variety of scenarios involving hatchery 
and natural production, ranging from natural production only, to production supplemented 
with hatchery adults or juveniles. The model is written in the SAS® programming language, 
which allows it to be used on a variety of computing systems and enhances flexibility in 
/ enstvzing model et 
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